JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOL10.1029/, 


Shocklets, SLAMS, and field-aligned ion beams in the terrestrial 
foreshock 
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Kasper 3 B.A. Maruca 3 M. Pulupa 4 C.S. Salem 4 M. Wilber 4 

Abstract. We present Wind spacecraft observations of ion distributions showing field- 
aligned beams (FABs) and large-amplitude magnetic fluctuations composed of a series 
of shocklets and short large-amplitude magnetic structures (SLAMS). The FABs are found 
to have T„ ~ 80-850 eV, V b /V ew ~ 1.3-2.4, T^/T,^ ~ 1-8, and n 6 /n 0 ~ 0.2-11%. Sat- 
uration amplitudes for ion/ion resonant and non- resonant instabilities are too small to 
explain the observed SLAMS amplitudes. We show two examples where groups of SLAMS 
can act like a local quasi-perpendicular shock reflecting ions to produce the FABs, a sce- 
nario distinct from the more-common production at the quasi-perpendicular bow shock. 

The SLAMS exhibit a foot-like magnetic enhancement with a leading magnetosonic whistler 
train, consistent with previous observations. Strong ion and electron heating are observed 
within the series of shocklets and SLAMS with temperatures increasing by factors >5 
and >3, respectively. Both the core and halo electron components show strong perpen- 


dicular heating inside the feature. [Date: 10/09/2012] 


1. Introduction 

Collisionless shock waves are a ubiquitous phenomena in 
space plasmas and are known to produce several different 
populations of reflected ion species, including: (1) field- 

aligned beams (FABs) [e.g. Greenstadt, 1976; Bonifazi and 
Moreno, 1981a, b], (2) intermediate ions [e.g. Paschmann 
et al, 1979], (3) diffuse ions [e.g. Gosling et al, 1978], (4) 
gyrating ions [e.g. Meziane et al., 1997], and (5) gyrophase- 
bunched ions [e.g. Gurgiolo et al., 1981]. The first three 
types are nearly gyrotropic and only distinguished by their 
pitch-angle distributions [e.g. Paschmann et al., 1981]. The 
difference between gyrating and gyrophase-bunched ions has 
been defined by their gyrotropy [e.g. Fuselier et al, 1986] 
and relative distance from the terrestrial bow shock [e.g. 
Meziane et al., 2001]. These ion distributions have been ex- 
amined in detail [e.g. Paschmann et al, 1981; Bale et al., 
2005; Eastwood et al., 2005, and references therein \ and all 
the ion species have been observed to be spatially well sep- 
arated [e.g. Meziane et al, 2011]. We will focus on FABs of 
finite spatial extent herein. 

The properties of reflected ion distributions depend 
strongly upon the shock geometry where the ion reflected 
from the shock surface. The shock normal angle, d Bn , is de- 
fined by the angle between the the average upstream mag- 
netic field and the local shock normal vector. The shock 
geometry is called quasi-perpendicular(-parallel) when d Bn 
>(<) 45°. Early observations [Greenstadt, 1976] suggested 
that FABs had their origin on field lines connected to a 
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quasi-perpendicular portion of the bow shock. Later studies 
found consistent results, namely that FABs are typically ob- 
served upstream of the quasi-perpendicular bow shock [e.g. 
Fuselier, 1995]. However, theory suggests that specular re- 
flection under quasi-parallel geometry results in a guiding 
center velocity directed upstream [e.g. Gosling et al., 1982], 
which has been supported by observations [e.g. Thomsen 
et al., 1990]. 

FABs are primarily composed of protons streaming along 
the ambient magnetic field away from the bow shock with 
temperatures (T 6 ) ~80-600 eV, densities (n t ) 1-10% of the 
ambient solar wind density (n 0 ), and beam speeds up to 
5 times the solar wind speed (V sm ) [Bonifazi and Moreno, 
1981a, b; Paschmann et al., 1981]. FABs often have strong 
temperature anisotropies with Tj_ it /T|| 6 > 4-9 [Paschmann 
et al., 1981]. They are typically observed in the absence of 
or near large magnetic fluctuations, but not simultaneously 
with these waves [Kis et al., 2007; Meziane et al., 2011]. 
Note that Meziane et al. [2004] observed FABs simultane- 
ously with gyrating ions near the ultra-low frequency (ULF) 
wave boundary, due to finite gyroradius effects. 

The terrestrial foreshock on the quasi-parallel side has 
a broad spectrum of large amplitude waves. This spec- 
trum consists of transverse Alfvenic waves, right-hand po- 
larized (in plasma frame) ULF waves near the ion gyrofre- 
quency (f ci ), compressional magnetosonic waves [Hoppe and 
Russell, 1983], magnetosonic- whistler mode waves [Hoppe 
et al., 1981], and an ensemble of higher frequency (f > 5-10 
Hz) waves up to the electron plasma frequency (f pe ) [e.g. 
Briand, 2009, and references therein]. We will focus on two 
specific types of waves in the foreshock, shocklets [Hoppe 
et al, 1981] and short large-amplitude magnetic structures 
(SLAMS) [Schwartz et al, 1992]. Both types are magne- 
tosonic in nature [Hoppe et al, 1981; Mann et al., 1994; 
Behlke et al, 2003; Hellinger and Mangeney, 1999], which 
means magnetic fluctuations, JB, are in phase with density 
fluctuations, <5n. The magnetosonic nature of shocklets and 
SLAMS causes them to dispersively radiate higher frequency 
electromagnetic whistler precursor waves as they steepen 
and they are always observed simultaneously with diffuse 
ion distributions [e.g. Wilson III et al., 2009, and refer- 
ences therein] . Hybrid simulations found that SLAMS could 
be produced by electromagnetic ion/ion beam instabilities 
[Onsager et al., 1991; Akimoto et al., 1993; Dubouloz and 
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Scholer, 1995; Hellinger and Mangeney, 1999]. Particle-In- 
Cell (PIC) simulations, however, have found that shocklets 
and SLAMS can result from a nonlinear interaction between 
gradients in the diffuse ion densities and the ULF wave field 
[e.g. Scholer et al., 2003; Tsubouchi and Lembege, 2004]. 
The exact nature of SLAMS formation is still unknown but 
it clearly requires nonlinear instabilities. 

SLAMS have been found to have wave vectors directed 
upstream away from the bow shock. Observations have 
shown that their phase speed is less than the solar wind 
speed but ^1-6 V A [e.g. Mann et al, 1994], where V A 
is the Alfven speed. Multi-spacecraft studies have found 
their convection speed (relative to the bow shock) to de- 
crease with increasing amplitude [e.g. Schwartz et al., 1992; 
Mann et al, 1994] and decreasing distance to the bow shock 
[e.g. Mann et al., 1994; Blanco-Cano, 2010]. This implies 
that larger amplitude SLAMS can stand against the inci- 
dent solar wind flow. Simulations have shown their wave 
vectors become more aligned with the local shock normal 
with decreasing distance to the shock [e.g. Dubouloz and 
Scholer, 1995], which supports the predictions of Schwartz 
and Burgess [1991]. Lucek et al. [2002] found that shocklets 
and ULF waves have larger scale sizes than SLAMS. Later, 
Lucek et al. [2004, 2008] determined that SLAMS scale sizes 
parallel to the shock normal were >1000 km and ~1300 
km parallel to the shock surface. They also found that the 
gradient scale length of the magnetic field was ~100 km. 
SLAMS and shocklets are more easily differentiated by their 
relative amplitudes. Shocklets exhibit weak magnetic com- 
pression with <5B/B 0 < 2, while SLAMS are much stronger 
with <5B/B 0 > 2 and often exceeding a factor of 4. 

Previous studies have found FABs near SLAMS, but in- 
strumental limitations prevented an examination of the evo- 
lution of the ion distributions across the SLAMS [Schwartz 
et al., 1992; Wilkinson et al., 1993]. Wilkinson et al. [1993] 
concluded that the SLAMS were generated by a right-hand 
resonant ion/ion instability [e.g. Onsager et al., 1991; Aki- 
moto et al., 1993]. They also concluded that the observed 
ion beams were most likely reflected off of a locally non- 
planar surface, like the SLAMS, but they do not exclude 
the bow shock as a source. Giacalone et al. [1993] found 
an increase in suprathermal ion pressure associated with 
SLAMS convecting back toward the bow shock. The pres- 
sure increase was observed on the side of the SLAMS facing 
the bow shock, suggesting the structures were preventing 
back-streaming suprathermal ions from escaping farther up- 
stream. These results coupled with observations showing 
gradient scale lengths ~100 km support the hypothesis by 
Mann et al. [1994] that SLAMS can act as efficient particle 
mirrors. 

In this paper we report the first high time resolution ob- 
servations of the evolution of FABs through large amplitude 
magnetic fluctuations, identified as shocklets and SLAMS, 
in the terrestrial foreshock. The ion beams are more intense 
on the upstream (sunward) side of the SLAMS, suggesting 
a local source. This is the first reported occurrence of a 
secondary mechanism for FAB production. The paper is 
organized as follows: Section 2 discusses the data sets and 
analysis techniques, Section 3.1 provides an overview of rel- 
evant parameters, Section 3.2 summarizes the particle dis- 
tribution observations, Section 3.3 justifies our conclusion 
that the group of SLAMS produce the FABs, and Section 4 
presents our conclusions. 

2. Data Sets and Analysis 

The magnetic field was obtained from the Wind dual, 
triaxial fluxgate magnetometers [Lepping et al., 1995] sam- 
pled at ~11 samples/s. Full 4-7T steradian low energy (<30 


keV) ion and electron distributions were obtained from the 
Wind/3DP EES A and PESA particle detectors [Lin et al, 
1995]. For more details about the analysis of data from 
the 3DP instrument, see Wilson III et al. [2010]. The solar 
wind velocity (V„„), ion thermal speed (V Ti , where V Ts 2 
= 2k B T s /m s for species s), and average ion temperature 
(Ti) were determined with the 3DP PESA Low and SWE 
Faraday Cups (FCs) [Ogilvie et al, 1995]. Absolute elec- 
tron densities were determined from the plasma line in the 
WAVES thermal noise receiver (TNR) instrument Bougeret 
et al. [1995] and used as an estimate of the solar wind ion 
density (n D ). 

To compare our observations to theory for the growth 
of SLAMS through a right-hand resonant ion/ion beam in- 
stability [e.g. Akimoto et al, 1993], we fit the FABs to bi- 
Maxwellians to determine estimates of the density (n 6 ), the 
temperature (T 6 ) or thermal speed (V Tb ), and beam speed 
(V 6 ). The maximum growth rate ("/max) for this instability 
is proportional to (n 6 /n 0 ) 1 ' /3 £2 ci , where Ll ci is the ion cy- 
clotron frequency. The corresponding wave number (k max ) 
normalized by the ratio "V A / LI ci , is proportional to V Ti /V;,, 
where V A is the Alfven speed. The corresponding real part 
of the frequency (cu rmox ) is proportional to (V-ri/V^ - 1) Ll ci , 
which can be related to the period of the SLAMS. The FABs 
are all observed when V b /V A 1- Under these conditions, 
the saturation amplitude of the right-hand resonant (and 
non-resonant) ion/ion beam instability is given by <5B/B„ 
~ (n i ,/2n 0 ) 1 / 2 (V !) /V A ) [e.g. Gary, 1991]. We will use these 
relationships to compare observations with theory. 

The wave vector, k, and the polarization with respect to 
the quasi-static magnetic field were determined using Min- 
imum Variance Analysis (MVA) [ Khrabrov and Sonnerup, 
1998]. The details of this technique are discussed in Wil- 
son III et al. [2009]. The uncertainties of this analysis are 
described by Khrabrov and Sonnerup [1998]. We calculated 
the angles between the wave vector and the local magnetic 
field ( d kB ) and the solar wind velocity ( 6 kv ). This analysis 
was only applied to the shocklets, SLAMS, whistler precur- 
sors, and foreshock fluctuations. 

We used two methods to determine the bow shock normal 
vector. The first involved the use of the Rankine-Hugoniot 
conservation relations [e.g. Koval and Szabo, 2008] with the 
parameters observed at the last crossing of the bow shock. 
From this we derive a single shock normal vector. The sec- 
ond method involved projecting the local smoothed mag- 
netic held vector onto the surface of a model bow shock 
[Slavin and Holzer, 1981]. Once we determined the shock 
normal vector, we were able to calculate the shock normal 
angle, d Bn . Due to the known errors inherent in each model 
and the dynamic nature of the bow shock, we use each shock 
normal estimate as a qualitative reference to assist in our 
analysis. 

3. Observations 

3.1. Foreshock Observation Overview 

Figure 1 shows the trajectory of the Wind spacecraft 
in the XY-GSE plane for the two dates examined. The 
projection of the average magnetic held vector in the fore- 
shock (green arrows) was determined by averaging the mag- 
netic held for the time period between the foreshock struc- 
tures (blue asterisks) and the foreshock boundary (green 
square). The projection of the average interplanetary mag- 
netic held (IMF, red arrows) was determined by averaging 
the smoothed magnetic held for the time period upstream 
of the foreshock boundary, which was 20:12-23:59 UT for 
2000-04-10 and 14:10-17:10 UT for 2002-08-10. We also in- 
cluded the projection of the model magnetopause and bow 
shock locations, where we adjusted the model parameters to 
match the observed crossings along the shown trajectories. 

The foreshock structures, identified as groups of shocklets 
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and SLAMS, were located at a GSE position of ~<+13.6, 
-12.4, +0.02> R e for the 2000-04-10 event (left panel) and 
^<+13.8, -3.0, +0.4> R E for the 2002-08-10 event (right 
panel). These positions correspond to distances of ~1.4 R E 
(^60 ion inertial lengths, c/oj pi ) and ~2.4 R E (^220 c/uj pi ) 
from the last bow shock crossings, respectively. Therefore, 
we needed to determine whether these structures were sim- 
ply due to an expansion of or close encounter with the ter- 
restrial bow shock. 

To exclude the possibility of a bow shock expansion, we 
examined data from the ACE, GOES 8 and 10, and Inter- 
ball spacecraft (courtesy of CDAWeb). We found no tran- 
sient sudden reduction/enhancement in solar wind pressure 
for either event and only small changes in the IMF orien- 
tation. Analysis of ACE and OMNI data showed a slight 
rotation in the Y-GSM IMF component for the 2002-08-10 
event. However, we do not believe it significant enough to 
explain a sudden expansion of >2 R E near the bow shock 
nose. Analysis of the GOES data showed no evidence of 
a sudden change in pressure applied to the magnetosphere, 
thus no change in magnetopause/bow shock location. Given 
that these structures are observed > 60 c/io pi upstream of 
the last bow shock crossing, we also do not believe these 
features to be features of a reformation process \e.g. Winske 
et al., 1990]. Therefore, we conclude that these structures 
are features of the terrestrial foreshock and not the bow 
shock. 

Next, we needed to determine the local geometry of the 
bow shock. Figure 2 plots an overview of the magnetic held 
measurements observed for the two groups of shocklets and 
SLAMS. We will focus on the 2002-08-10 event herein. The 
blue line in Figures 2E and 2J represents the d Bn calcu- 
lated from our first method and the magenta line was cal- 
culated using the second method, both discussed in Section 
2. Note that the second method used a smoothed average 
magnetic held, instead of the HTR MFI data, and gaps indi- 
cate regions not magnetically connected to the model shock 
surface. The Rankine-Hugoniot solutions gave us shock nor- 
mal vectors of ~<+0.750, -0.621, -0.134> for the 2000-04- 
10 event and ~<+0.987, -0.158, +0.002> for the 2002-08-10 
event. We estimated the upstream foreshock average mag- 
netic Held, B„, using the Wind/MFI observations between 
14:20-15:22 UT for the 2000-04-10 event and between 11:21- 
12:49 UT for the 2002-08-10 event. The average GSE vec- 
tors were found to be B 0 ~ (+3.60,-2.99,-1.47) nT and ~ 
(—4. 32, +0.65, —0.41) nT, respectively. These estimates with 
the above shock normal vectors give 0 Bn ~ 14° and ~ 6°, 
respectively. Therefore, the Wind spacecraft was immersed 
within the quasi-parallel region of the terrestrial foreshock 
for the time of interest. 

The foreshock structures, marked by vertical red lines 
in Figures 2A-B and 2F-G, were composed of a series of 
compressive magnetosonic waves (<5B in phase with density, 
<5n, fluctuations) identified as shocklets and SLAMS shown 
in Figures 2C-D and 2H-I, respectively. The SLAMS were 
observed to have: (1) mixtures of right- and left-hand polar- 
izations (spacecraft frame); (2) very oblique (d kB > 55° and 
e kV > 40°) propagation; and (3) <5B/B„ > 2-6, consistent 
with previous observations [Schwartz et al, 1992; Wilkin- 
son et al., 1993]. The 2002-08-10 event (Figure 2H-I) shows 
an isolated SLAMS near 12:52:20 UT. The 2000-04-10 event 
(Figure 2C-D) did not show a similar structure. The impor- 
tance of this difference will be discussed in the next section. 

Both groups of SLAMS have higher frequency fluctua- 
tions on their leading/upstream (i.e. the right-hand side 
of Figures 2C-D and 2H-I) edges, consistent with whistler 
mode waves. The characteristics of the waves immediately 
upstream(sunward) of the the steepened edges are consis- 
tent with previous observations of whistler precursors [ e.g . 
Wilson III et al., 2009]. The whistler amplitudes and beam 
intensity decrease away from the leading edge of the group 


of SLAMS and eventually the whistlers disappear when the 
FABs disappear. However, we cannot dehnitively show that 
the two phenomena are causally related because whistler 
modes have been observed in the absence of FABs. While 
simulations have found that rehected ions can provide free 
energy to enhance already present whistler precursors [e.g. 
Scholer et al., 2003], supported by recent observations [ Wil- 
son III et al., 2012], their primary source is thought to be 
dispersive radiation [e.g. Sundkvist et al, 2012]. The source 
of the whistler mode waves are beyond the scope of this 
paper and we will not discuss them further. 

3.2. Particle Distributions 

We examined the effects caused by the series of shocklets 
and SLAMS on the ion and electron distribution functions 
for the two foreshock passes. The group of SLAMS created 
a rarefaction region behind the structures with a strong de- 
flection of the solar wind core, analogous to the wake created 
by an obstacle in a huid how. The SLAMS caused strong 
anisotropic heating in the low energy (<1.1 keV) electrons 
and ions (<10 keV). 

Figure 3 shows an example ion distribution (in the bulk 
How frame) plotted as contours of constant phase space den- 
sity versus velocity (axes range from ±1500 km/s) contain- 
ing a FAB. To aid the reader, we have inserted ellipses to 
show the location of the FAB (black) and solar wind beam 
(purple). In the top distribution, one can see the FAB (near 
~400 km/s) is moving anti-parallel to B 0 , which corresponds 
to the sunward direction for this event. Note that lowest 
energy bin for PESA High for both events was ^80 eV (or 
~125 km/s proton). This figure is used for illustrative pur- 
poses and the format for all distributions shown in Figure 4 
are consistent with the top contour plot in Figure 3. This 
distribution corresponds to panel B in Figure 4. 

Figure 4 shows the HTR MFI data and select PESA High 
distribution functions for the time range corresponding to 
Figure 2H-I. The solar wind core is clearly identified in the 
center of sequential panels S to AE and separated from the 
FABs seen near ^500-900 km/s in panels S-AB and ^400 
km/s in panel B (see Figure 3). The FAB in panel B is at 
lower speeds and comparison with foreshock distributions 
(e.g. panel A) improves contrast for better identihcation. 
The FABs on the downstream (earthward) side of the group 
of SLAMS are weaker for both events. Hot diffuse ions (i.e. 
nonthermal tail observed above ~800 km/s) are observed 
continuously between 12:50:13 UT and 12:51:45 UT, simul- 
taneous with the SLAMS. PESA Low and SWE distribu- 
tions showed that the core ions experienced very strong heat- 
ing in this region, strong enough to be observed by PESA 
High in panels D-R in Figure 4. The observed durations 
(At FAB ) of FABs and effects on core particle populations 
are consistent with previous observations [ Schwartz et al., 
1992; Wilkinson et al., 1993; Meziane et al., 2004], 

The FAB properties, summarized in the top half of Ta- 
ble 1, show the range of values found from our ht results. 
The values for the 2002-08-10 properties are for FABs ob- 
served between 12:51:51-12:52:48 UT and for FABs satisfy- 
ing V b /V A >4 (only one FAB had V b /V A < 4). These prop- 
erties are consistent with previous observations upstream of 
the quasi-perpendicular bow shock [e.g. Bale et al, 2005, 
and references therein}. The bottom half of Table 1 shows 
the instability analysis results. For reference, the magnetic 
held values used to produce the green arrows in Figure 1 
correspond to Q ci ~ 0.44 radians/s for the 2000-04-10 event 
and Q ci ~ 0.52 radians/s for the 2002-08-10 event. 

Note that FABs were observed both downstream (earthward) 
and upstream (sunward) of the group of SLAMS for both 
events (2000-04-10 event not shown). The FAB intensity 
was found to be greater on the upstream(sunward) side of 
the group of SLAMS but weaker on the upstream(sunward) 
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side of the isolated SLAMS (^12:52:20 UT) shown in Figures 
2H-I. Figure 5 shows a zoomed in view of the 2002-08-10 
event where FABs are observed upstream(sunward) of the 
group of SLAMS. The timestamps for the data points were 
determined from the average of the start and end times of 
each PESA High distribution. One can see that n 6 and n 6 /n„ 
decrease until ~12:52:15 UT then rise from ^12:52:19 UT 
(Figure 4W) to ~12:52:22 UT (Figure 4X). Note that the 
local peak at ~12:52:22 UT corresponds to the ramp of the 
isolated SLAMS and that n b and n 6 /n 0 decrease thereafter. 
The bottom two panels in Figure 5 show that V b /V 3 v, and 
Vi ,/Y a increase steadily over this same interval. 

The decrease in n b /n 0 with increasing distance up- 
stream(sunward) from the group of SLAMS (shown in Fig- 
ures 2H-I) suggests a local source. The large scale size 
(~4000 km) and comparable gradient scale size (~100 km) 
of SLAMS to the average ion inertial length (~100 km) al- 
lows them to act as efficient magnetic mirrors when grouped 
together. Isolated SLAMS, however, are known to have 
lower amplitudes and thus lower phase speeds [ e.g . Schwartz 
et al., 1992], which results in a smaller velocity between inci- 
dent solar wind ions and the SLAMS. The increase in n 6 /n 0 
on the downstream side of the isolated SLAMS suggests 
that this structure is preventing some beam ions from es- 
caping further upstream, consistent with previous observa- 
tions [ Giacalone et al, 1993]. The local peak in n b /n 0 in the 
upstream(sunward) edge of the isolated SLAMS shows that 
it can reflect particles as well, but not as effectively as the 
group of SLAMS. The reasons for why an isolated SLAMS 
cannot reflect as well as the group of SLAMS may be: (1) 
isolated SLAMS are too narrow spatially; (2) they have 
smaller <5B/B 0 ; (3) they do not propagate as fast causing 
their phase velocity relative to the bulk flow to be smaller; 
and/or (4) less rotation and turbulence in the magnetic field. 

If the local magnetic field cannot connect the spacecraft 
with the bow shock without going through the group of 
SLAMS, then one would not expect to observe the coherent 
FABs shown in Figure 4. Therefore, if the group of SLAMS 
are between the spacecraft and the bow shock along an aver- 
aged magnetic field direction, then they must be the source 
of the FABs. 

3.3. FAB Source 

In this section we will explain why we believe the FABs 
originate from the SLAMS and not the bow shock. First 
we need to eliminate the possibility that the FABs produce 
the SLAMS by comparing observations to theory for ion/ion 
beam instabilities. Then we will use geometry to argue that 
the group of SLAMS create a magnetic barrier between the 
spacecraft and the bow shock. Lastly we will summarize our 
results and present our conclusions. 

The observations show (Figures 2C-D and 2H-I) that 
the duration of individual SLAMS are longer for the 2002- 
08-10 event than the 2000-04-10 event. The predicted wave 
periods (w rma x _1 ) from Table 1 were ~1.6-3.7s for the 2000- 
04-10 event and ~5. 4-25. Is for the 2002-08-10 event. The 
correspondence between the predicted wave periods and ob- 
served durations might initially suggest that the FABs are 
causing an ion/ion beam instability resulting in the ob- 
served SLAMS. However, this correspondence requires that 
the structures convect with the bulk flow and that they do 
not propagate in the plasma frame. 

The average upstream (sunward) V a „ (^460 km/s) is 
comparable for both events and the distribution of 6 kv are 
roughly the same. Therefore, we expect comparable Doppler 
effects for both events. Previous observations have shown 
that each individual SLAMS has a normal scale length L n 
~ 1000 km and transverse scale length L, ~ 1300 km [Lucek 
et al, 2008]. These scales correspond to convection time 
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scales of L„/V sm ~ 2.2 seconds per SLAMS and L t /V s „ ~ 
2.8 seconds per SLAMS. These values do not seem to cor- 
relate with the observed wave periods, but we have ignored 
the phase speed of the SLAMS in these estimates. 

Recall that the phase speed of SLAMS was determined to 
be V pfe ~l-6 Va [Mann et al, 1994] and they are known to 
propagate nearly along the bow shock normal vector [Lucek 
et al., 2008]. The geometry for both events is such that the 
spacecraft frame speed of the SLAMS is V 0 = | V s „ + ~V ph 
| < V.„. For example, if we use V pfe = 4 V A along the shock 
normal, then V„ ~ 200-250 km/s for the 2000-04-10 event 
and V 0 ~ 80-230 km/s for the 2002-08-10 event. These new 
speeds change the convection time scales to L„/V 0 ~ 4.1- 
5.0 seconds per SLAMS for the 2000-04-10 event and L„/V 0 
~ 4.4-12.6 seconds per SLAMS for the 2002-08-10 event. 
Therefore, the longer observed durations for the SLAMS for 
the 2002-08-10 event can be explained by a slower convec- 
tion speed of the structures. 

Now we examine whether the beams could produce wave 
amplitudes large enough to match the observations. When 
we examined the parameters in Table 1 we used the results 
to estimate the saturation amplitude of the right-hand reso- 
nant and non-resonant ion/ion beam instabilities for V b /V A 
1 [e.g. Gary, 1991]. In this limit we find 5B/B 0 ~ 0.4 ± 
0.1 for the 2000-04-10 event and <5B/B„ ~ 1.3 ± 0.1 for the 
2002-08-10 event. These values represent the mean plus or 
minus the standard deviation of the mean for all the beam 
distribution fit results for each event. Examination of the 
SLAMS in Figures 2C-D and 2H-I show (5B/B„ ~ 3-6, much 
larger than theory predicts. More importantly, the FABs are 
not observed inside the group of SLAMS. These results ar- 
gue against the ion/ion beam instability as a source for the 
SLAMS. 

Figure 6 is an illustrative cartoon that we will use to ar- 
gue that the group of SLAMS, referred to as the obstacle 
for brevity, in each event is the source of the FABs. If we 
assume the obstacle is being converted with the solar wind 
at roughly V.,„, then the spacecraft will be effectively sta- 
tionary. Therefore, the path of the spacecraft through the 
obstacle is anti-parallel to V 3U) with a length of L s = V sw 
At sc . The amount of turbulence and rotation in the mag- 
netic field observed through the obstacle should serve as a 
magnetic barrier from particles leaking into the upstream 
from behind the obstacle. Therefore, the spacecraft will be 
in the “magnetic shadow” of the obstacle, with respect to 
the terrestrial bow shock, for a distance L shadow . 

Recall that each individual SLAMS has scale lengths ~ 
1000 km. It takes At ac ~ 136s(108s) to traverse the ob- 
stacle and we observe FABs for At F As ~ 10s(57s) for the 
2000-04-10(2002-08-10) event. The duration the spacecraft 
spends within the obstacles (At sc ) suggests b<Ca. The av- 
erage complementary angle between V sm and B 0 is < 25° for 
both events. Therefore, any trajectory through the obstacle 
above the horizontal in Figure 6 will result in the spacecraft 
being in the “magnetic shadow” of the obstacle for at least 
At sc after exiting on the upstream (sunward) side. 

If the spacecraft trajectory is below the horizontal in Fig- 
ure 6 then we expect the duration in the “magnetic shadow” 
of the obstacle to be shorter than for trajectories above the 
horizontal. However, even for the trajectories resulting in 
significantly shorter durations of “magnetic shadow,” the 
observations of FABs in such close proximity to the obstacle 
argues against a bow shock source. 

4. Discussion and Conclusions 

This study presents observations of field-aligned ion 
beams (FABs) near large amplitude magnetic held fluctu- 
ations in the terrestrial foreshock by the Wind spacecraft. 
We examined two foreshock events, on 2000-04-10 and 2002- 
08-10, composed of groups of shocklets and short large- 
amplitude magnetic structures (SLAMS). The 2002-08-10 
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event exhibited an isolated SLAMS upstream(sunward) of 
the group which allowed us to examine the evolution of the 
FAB through the structure. 

The FABs propagate (in the plasma frame) away from 
the bow shock toward the upstream(sunward) side of the 
SLAMS. They had T b - 80-850 eV, V b /V ew ~ 1.3-2.4, 
Tx.t/Ty^ ~ 1-8, and n b /n a ~ 0.2-11%, consistent with pre- 
vious observations [e.g. Bale et al, 2005, and references 
therein ]. While ion beams have been previously observed 
near SLAMS [Schwartz et al, 1992; Wilkinson et al., 1993], 
no previous reports have shown the evolution of ion distri- 
butions through SLAMS and showed the SLAMS to be the 
source of FABs. 

We observed peak values of n b /n„ near the up- 
stream(sunward) edge of the group of SLAMS in both 
events. Another peak was observed on the immediate up- 
stream(sunward) edge of the isolated SLAMS near 12:52:20 
UT in the 2002-08-10 event. Just downstream (earthward) 
of the isolated SLAMS n b /n a increased suggesting the 
SLAMS was preventing beam ions from escaping up- 
stream(sunward). The FABs for this event also had larger 
values of n 6 , T 6 , T Xii ,/T|| ib , and V b /V A - These differences 
are probably due to larger amplitude of the SLAMS in the 
2002-08-10 event. 

The source of the whistler precursors is beyond the scope 
of this manuscript. The SLAMS, however, are thought to 
be driven by gradients in the diffuse ion density [e.g. Sc- 
holer et al, 2003]. We observed diffuse ions throughout each 
group of SLAMS and we did not observe FABs within either. 
Therefore, we believe the diffuse ions, not the FABs, are re- 
sponsible for the SLAMS. 

In summary, we argue that the SLAMS are the source 
of these FABs for the following reasons: (1) the decrease in 
n b /n 0 with increasing distance from the upstream(sunward) 
side of the group of SLAMS; (2) the evolution of the FABs 
across the isolated SLAMS shows a local peak in n 6 /n 0 
on the immediate upstream(sunward) edge; (3) FABs are 
not observed within the group of SLAMS and their peak 
intensity is at the immediate upstream (sunward) edge of 
the group of SLAMS; (4) the predicted saturation ampli- 
tudes for the ion/ion beam instability are much smaller 
than observed; (5) observed differences in duration can be 
explained by differences in convection speed of the struc- 
tures; and (6) geometry shows that the spacecraft is in the 
“magnetic shadow” of the group of SLAMS at the immedi- 
ate upstream(sunward) edge. Therefore, at least the FABs 
observed immediately upstream(sunward) of the group of 
SLAMS cannot have a bow shock source and our results ar- 
gue that it is likely the rest do not originate from the bow 
shock either. 

The similarity of the field-aligned ion beams with those 
observed upstream of the quasi-perpendicular shock sug- 
gest the groups of SLAMS are acting like a local quasi- 
perpendicular shock. This is not to say that all field-aligned 
ion beams are produced by SLAMS. However, it is possi- 
ble the beams are produced in a manner similar to those 
found in simulation results, showing a positive feedback loop 
between the waves and reflected particles [e.g. Scholer and 
Burgess, 1992]. 

In conclusion, we show the first direct evidence that 
groups of SLAMS can locally produce field-aligned ion 
beams. 
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Table 1. Instability Analysis Results 


Range of Beam Parameters 

Date 

T, (eV) 

Tj_,!,/T|| ib 

n 6 /n 0 (%) 

V„/V 3 „ 

v b /v A 

v„/v Tb 

2000-04-10 

80-170 

1. 1-4.4 

0.2-1. 6 

1. 5-2.0 

7. 5-8. 5 

3. 7-5.0 

2002-08-10 

175-850 

2. 3-8. 2 

0.3-11.0 

1.3-2. 4 

9.1-18.1 

2. 1-3.7 

Instability Results 

Date 

At fab ( s ) 

'Y max / ^ ci 

k max ^ A / ^ci 

bJrmax /^ci 

'Y max (s) 

2000-04-10 

-10 

0.10-0.20 

0.18-0.32 

0.34-1.43 

9.10-19.6 

2002-08-10 

-57 

0.11-0.38 

0.064-0.12 

0.043-0.23 

2.62-11.3 



Figure 1. The trajectory (magenta line) of the Wind 
spacecraft in the XY-GSE plane for the 2000-04-10 (left 
panel) and 2002-08-10 (right panel) events. The trian- 
gle corresponds to the start time shown in the plot ti- 
tles and the diamond the end time. The solid blue arc 
shows the projection of the model magnetopause location 
[ Kawano et al, 1999] and the solid black arc the model 
bow shock location [Slavin and Holzer , 1981]. The green 
arrows show the projection of the average magnetic held 
vector in the foreshock and the red arrows show the corre- 
sponding projection for the solar wind. The green square 
marks the approximate foreshock boundary location. 
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Figure 2. The top half of the figure shows three second 
resolution of the magnitude and the normal incidence 
frame (NIF) Sundkvist et al. [e.g. 2012] components of 
magnetic field data from the Wind spacecraft on 2000-04- 
10 (A-B) and 2002-08-10 (F-G) each with three vertical 
lines that indicate the magnetopause crossing (green), the 
last bow shock crossing (blue), the foreshock boundary 
(purple), and the red lines show the time periods for pan- 
els C-E and H-J. The tick mark labels at the bottom of 
these two panels are: UT time, and the Wind spacecraft 
radial distance (Re), GSE longitude (degrees), and GSE 
latitude (degrees). Every panel has the same format, but 
the bottom two panels show the HTR MFI data and the 
shock normal angle for model (magenta) and Rankine- 
Hugoniot (blue) solutions. 


RH Solution Fluxgate Magnetometer Fluxgate Magnetometer 

Model Solution ~11 Samples/s ~0 .33 Samples/s 
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Wind 3DP PESA High Burst 
2002-08-10/12:50:04 - 12:50:07 
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Figure 3. An example PESA High Burst distribution 
observed at 2002-08-10/12:50:04 UT. The contour plots 
show contours of constant phase space density (uniformly 
scaled from lxlO -13 to lxl0~ 9 s 3 cm -3 km -3 , where red 
is high) projected onto three different planes defined by 
the shaded region in the coordinate axes shown in right- 
hand column. The diagonal axis defines the normal to the 
plane shown for each corresponding contour plot. Pro- 
jected onto each contour are ~V aw (black arrow), the sun 
direction (blue arrow), and the Earth direction (magenta 
arrow). The FAB (black ellipses) and core (purple el- 
lipses) are labeled. This distribution corresponds to panel 
B in Figure 4. 
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Figure 4. Selected PESA High Burst distributions 
shown for the time range shown in Figure 2H-J. The 
ion distribution plots are the same format as the top left 
panel of Figure 3. Panels S-AE have circles of constant 
energy at 500, 700, 900, and 1100 km/s. Ion beams are 
clearly identified in panels B and S-AB. 
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Figure 5. This figure shows the FAB parameters versus 
time with the magnitude of the magnetic held for the 
2002-08-10 event. The panels, from top-to-bottom, show: 
magnetic Held magnitude, B„ (nT); FAB number density, 
n b (cm -3 ); ratio of FAB to total number density, n b /n 0 ; 
ratio of FAB drift speed to bulk How speed, V o6 /V au ,; and 
ratio of FAB drift speed to local Alfven speed, V 0 i ,/V a- 



Figure 6. A schematic cartoon used to illustrate how 
the series of SLAMS in Figures 2C-E and 2H-J can block 
the spacecraft (SC) from “seeing” the bow shock along 
magnetic field lines. In this example, the sun is to the 
right and Earth to the left. 


